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ABSTRACT 

Extended nebular filaments are seen at optical wavelengths in NGC 1275, the 
central galaxy in the Perseus cluster. The agents responsible for the excitation of 
these filaments remain poorly understood. In this paper we investigate possible 
mechanisms for powering the filaments, using measurements from an extensive 
spectroscopic data set acquired at the Lick Observatory 3-m Shane telescope. 
The results show that the filaments are in an extremely low ionization and exci- 
tation state. The high signal-to-noise ratio of the spectra allows us to measure 
or place sensitive upper limits on weak but important diagnostic lines. We com- 
pare the observed line intensity ratios to the predictions of various ionization 
models, including photoionization by an active galactic nucleus, shock heating, 
stellar photoionization, and photoionization by the intracluster medium. We also 
investigate possible roles for cluster extreme-ultraviolet emission, and filtering of 
cluster soft X-ray emission by an ionized screen, in the energetics of the filaments. 
None of these mechanisms provides an entirely satisfactory explanation for the 
physical state of the nebulae. Heating and ionization by reconnection of the in- 
tracluster magnetic field remains a potentially viable alternative, which merits 
further investigation through Faraday rotation studies. 



Subject headings: galaxies: clusters: individual (Perseus) — cooling-flows — 
galaxies: individual (NGC 1275) — galaxies: ISM 



1. Introduction 



Extended nebular filaments are commonly seen in association with the central galaxies 
of galaxy clusters. Clusters that feature a massive central galaxy are also often identified 
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as having a cooling flow, based on cooling-timescale arguments or other evidence (Fabian, 
Nulsen, & Canizares 1984). This association has prompted spccTilation that the optical 
nebulosity is linked causally or energetically to the intracluster medium (ICM), and hence 
may be an important diagnostic of the cooling-flow phenomenon. As a result, several previous 
studies have investigated optical "cooling-flow fllamcnts" in attempts to understand their 
physical properties. These efforts have helped delineate the observed characteristics of the 
nebulosity, and have prompted a variety of theoretical explanations for what is seen. 

Optical spectroscopy reveals that the fllaments emit strongly in low-ionization forbidden 
lines, thus resembling Low Ionization Nuclear Emission-Line Regions (LlNERs; Hcckman 
1980). Heckman et al. (1989) also reported evidence that cooling-flow fllaments can be 
separated into two characteristic groups based on their line-intensity ratios. Typical ratios 
for lype I systems are [N II]A6583/Hq; ^ 2, [S II]A6717/Hq; fa 0.7, and [O I]A6300/Hq; ^ 0.2; 
these objects tend to have low X-ray and Ha luminosities and are associated with low 
cooling rates and small optical fllaments. In contrast, typical ratios for Type II objects 
are [N II]A6583/Hq; 1, [S II]A6717/Ha fa 0.4, and [O I]A6300/Ha ^ 0.3; these sources 
have high X-ray and Ha luminosities and are associated with high inferred cooling rates and 
extended optical fllaments. Theoretical efforts to interpret the optical spectra have invoked 
various combinations of cooling or irradiation from the ICM, shocks, active galactic nucleus 
(AGN) emission or outflows, magnetic reconnection, and hot stars. However, none of these 
scenarios has proven entirely satisfactory in explaining the observations; consequently, the 
origins and energetics of these systems remain poorly understood. 

An important example of optical fllaments in apparent association with a cooling flow 
is found in NGC 1275 (Kent & Sargent 1979), the central elliptical galaxy in the Perseus 
cluster (Abell 426). The nebular gas extends to at least 20 kpc from the center of the galaxy 
(assuming Ho=75 km Mpc~^), with a complicated spatial structure; the optical spec- 
trum fits within Heckman et al.'s Type II category. The surrounding intracluster medium is 
believed to have a cooling flow, with mass deposition rates of up to 130 Mq yr~^ (Fabian, 
Cowie, & Grindlcy 1981; Allen et al. 1999). NGC 1275 features additional unusual char- 
acteristics that may be relevant to interpreting the extended nebular emission. The central 
elliptical galaxy is seen at a recession velocity of ~ 5200 km s~^, but is partially obscured by 
a smaller foreground galaxy at a velocity of ~ 8200 km s^^ (Minkowski 1957); the redshift 
difference is consistent with the velocity dispersion observed for this cluster (e.g., Kent & 
Sargent 1983). NGC 1275 is also host to a strong AGN and associated radio source, 3C 84, 
which emits jets that may be interacting with the ambient gas (Pedlar, Booler, & Davies 
1983; Pedlar et al. 1990). 

In this paper, we report the results of a detailed study of the physical properties of the 
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optical filaments associated with NGC 1275. This work makes use of extensive spectroscopic 
observations to quantify two-dimensional aspects of the observed emission, and also places 
sensitive limits on weak spectral features of diagnostic interest. The results can be used to 
place stringent limits on the role of various agents proposed for heating or ionization of this 
plasma. 

2. Observations and Data Reduction 

Long-slit spectra of the filaments of NGC 1275 were acquired in 10 observing runs 
during 1988 - 1990, using the Shane 3-m telescope + UV Schmidt CCD spectrograph at Lick 
Observatory (Miller & Stone 1987). The resulting data cover the wavelength range 3396 - 
6940 A in three bandpasses, spanning 3396 - 5000 A, 4695 - 6280 A, and 6130 - 6940 A. 
The slit was 2.5" wide and 132" long, with spatial sampling of 0.66" pixel^^. Observations 
were acquired with the slit centered on the nucleus and oriented at four different position 
angles (PA = 62°, 77°, 283°, and 313°) selected to lie along filaments with the highest 
surface brightness in Ha (see Figure 1). Typical atmospheric seeing was 2"and not all 
observations were done under strictly photometric conditions. The spectral resolution was 
~3.5 A longward of 6130 A and ~7 A elsewhere. The giant H II region detected by Shields & 
Filippenko (1990a) lies at PA = 62°. Line emission from the high- velocity system {cz ^ 8200 
km s~^) appears primarily northwest of the nucleus (e.g., Unger et al. 1990); it is visible in 
our spectra acquired at PA = 313°, and to a lesser extent at PA = 283°. 

The data were reduced using standard techniques, as implemented in IRAF []. Spectra 
acquired at the same position angle and wavelength setting were averaged after flux and 
wavelength calibration and sky subtraction. Estimates of the sky contribution were derived 
from regions near the ends of the slit, where no nebular emission was evident and the con- 
tribution of galaxy starlight was minimal. Measurements of filament emission were obtained 
for regions where Ha emission was detected in a visual inspection of the 2-dimensional (2-D) 
spectra. 

Emission-line strengths for the stronger lines were measured along the slit on a pixel-by- 
pixel basis, using an automated algorithm. For each resulting 1-D spectrum, the continuum 
in each wavelength setting was approximated by a low-order polynomial fit to emission-free 
bandpasses. After subtraction of the continuum, candidate emission lines were measured 



^The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA), 
under cooperative agreement with the National Science Foundation. 
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by fitting Gaussian profiles to all features detected at > 3a significance. The lines were 
subsequently identified, and spurious features rejected, by comparison of the mcasiucd list 
with an input list of line wavelengths, shifted by an a priori estimate of the emission rcdshift; 
this redshift value was measured interactively for Ha in the spectrum of the filament coaddcd 
along the slit. Lines were identified if the predicted and observed wavelengths agreed to 
within 5 A. This bound was set to allow for the correct identification of lines in the presence 
of substantial (< 300 km s~^) velocity gradients that are observed in this gas (e.g., Rubin et 
al. 1978). The fine flux measured by this method is \/2T^a\Fp, where a\ is the line profile 
dispersion as determined by the Gaussian fit and is the peak height of the emission line 
above the continuum. 

Estimation of line-intensity ratios for the NGC 1275 nebulosity requires correction for 
reddening. The filaments are subject to a Galactic foreground extinction of Ay — 0.82 
mag (Lockman & Savage 1995) as well as internal extinction, some of which may be due to 
the high- velocity system. The hmited signal-to-noise ratio (S/N) of many of the pixel- by- 
pixel spectral extractions translates into significant uncertainties in a reddening value based 
on Balmer-line ratios, and as a result we derived an average reddening correction for each 
filament (i.e., region of nebular emission seen on one side of the nucleus at a given PA). We 
obtained total extinction values from the measured Ha/H/? line- intensity ratio, assuming an 
intrinsic value of 2.86 (Hummer & Storey 1987) and the extinction curve of Cardelli, Clayton, 
& Mathis (1989). The results were compared with extinction values implied by fine ratios 
employing H7 and H(5, and were found to be consistent, with the exception of measurements 
for PA = 313°, where the limited S/N made the high-order Balmer lines unmeasurable, and 
for PA = 283°. The line ratios at the latter PA show evidence of underlying stellar absorption, 
which is not surprising in light of earlier reports of early-type stellar contributions to the 
integrated light seen in some portions of NGC 1275 (e.g., Wirth ct al. 1983). The infiuence 
of the stellar features is apparently strongest at this PA since the emission equivalent widths 
are relatively weak. For this filament we determined a consistent reddening value for the 
higher Balmer transitions subject to absorption with constant equivalent width. The total 
extinctions (including the Galactic contribution, and neglecting the redshift difference of the 
absorbers) calculated by these means are Ay — 1.0±0.2, 0.8±0.2, 2.3±0.5, and 1.1±0.4 mag 
for filaments at PA = 62°, 77°, 283°, and 313°, respectively. The larger value of extinction 
towards PA = 283° may be caused by the high- velocity system (Unger et al. 1990). The 
B-band atmospheric feature appears slightly longward of redshifted [S II]A6731, but does not 
affect measurements of that line. 

In addition to studying the spatially resolved emission at each PA, wc generated a 
composite spectrum of the filament system in order to study weak emission features. For 
each PA, the 2-D spectra were averaged over intervals of high S/N filament emission (S/N 
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;^ 5 for Ha). For PA = 62°, this criterion translated into an extracted region 14" wide, 
centered 23" away from the nucleus. The PA = 77° and PA = 283° extractions were 18.5" 
wide, centered at 21" and 20" away from the nucleus, respectively. The extraction for PA 
= 313° spanned 33", centered at 26" from the nucleus; we ultimately decided not to include 
this extraction in the composite spectrum because the emission-line equivalent widths were 
generally less than for the other PAs, with a correspondingly lower S/N. The resulting 1- 
D spectra were shifted to the rest wavelength scale, as determined from the redshift of the 
Balmer lines, and the spectra for different PAs were then combined to a single overall average. 
The relative velocity shifts between the filament extractions were less than 90 km s~^. The 
final composite spectrum (Figure 2) represents a total integration time of 10 hours in each 
wavelength setting. 

In order to improve the accuracy of the continuum subtraction for the composite spec- 
trum, we experimented with modeling the stellar continuum with templates derived from 
observations of emission-free galaxies. A reasonable correspondence in terms of the strength 
of stellar features in NGC 1275 was found with the spectrum of NGC 205. The stellar velocity 
dispersion of NGC 1275 is 250 km s'^ (Nelson & Whittle 1995) while that of NGC 205 is 50 
km s~^ (Held, Mould, & de Zeeuw 1990). We consequently convolved the NGC 205 spectrum 
with a Gaussian having a = 245 km s~^ in order to match the dispersion seen in NGC 1275. 
To optimize the template match, a fit with minimization was carried out using SPECFIT 
as implemented in IRAF (Kriss 1994), varying the template redshift, scahng, and reddening. 
A best fit was obtained with Ay — 0.6 ±0.1 mag. The Galactic extinction for NGC 205 
is Ay — 0.10 mag (Burstein & Heiles 1984); if the two galaxies were perfectly matched in 
spectral properties, this would imply a total extinction for NGC 1275 of Ay = 0.7 ± 0.1 
mag, in agreement with the Galactic foreground estimate. This value is somewhat less than 
is obtained from Balmer emission lines for the filaments; the disparity may stem from real 
differences in the reddening of starlight and nebular emission, and/or differences in stellar 
populations for NGC 1275 and NGC 205. The spectral regions over which we performed the 
continuum fitting do not appear to be affected strongly by metallicity differences between 
NGC 1275 and NGC 205. 

For the stronger emission features, line strengths measured after subtraction of this 
continuum fit were consistent with the results obtained from a smooth fit to the continuum, 
as employed for the pixel-by-pixel extractions. Use of the galaxy template was important 
for obtaining accurate values or limits for the weaker features, however, and for some lines 
the template fit was locally optimized (by adjusting the scaling) in order to match the 
observed continuum. The line fluxes obtained after subtraction of the continuum template 
are listed in Table 1, and were obtained from Gaussian profile fits as described above. Error 
bars listed for the detected lines represent the root-mean-square uncertainty per pixel (cTc) 
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in the local continuum, summed in quadrature over a width of 6a x, i.e., y/6acCrx. For 
lines that were not detected, we list 3cr upper limits of \/27rax{3ac), where a\ is the width 
of either H/3 or Ha, depending on which of these Balmer lines is closer in wavelength. 
For the composite spectrum, the Ha/H/? ratio implies a total extinction of Ay = 1.21 
mag, yielding an intrinsic value of Ay = 0.39 mag. The higher-order Balmer lines in the 
continuum-subtracted spectrum showed evidence for being over-corrected for the underlying 
stellar absorption: {E.6/B.P) observed > (H(5/H/9)caseS- The implied disparity between the 
NGC 1275 and template spectra does not significantly influence the estimate of Ay obtained 
from Ha/up, however, due to the large equivalent widths of these emission features. 



3. General Results 

3.1. Line-Ratio Diagrams 

An informative overview of the emission properties of the NGC 1275 fllaments can 
be obtained from standard line-ratio diagrams employed for nebular classification (Baldwin, 
Phillips, & Terlevich 1981; Veilleux & Osterbrock 1987). Figure 3 shows three such diagrams 
along with measurements from the pixel-by-pixel extracted spectra for NGC 1275. The 
observations have been corrected for reddening, as described in §2, and the use of intensity 
ratios of lines close in wavelength makes these diagrams insensitive to any residual reddening 
uncertainties. Loci commonly adopted for classification of H II regions, Seyfert nuclei, and 
LINERs are also indicated (e.g.. Shields & Fihppenko 1990b). 

As found in previous work, the filaments in NGC 1275 are characterized by low-ionization 
emission, resembling that of LINERs (e.g., Heckman et al. 1989). An additional dia- 
gram illustrating the ionization state of the gas is shown in Figure 4, which is a plot of 
[O l]A6300/[O III]A5007 versus [O Il]A3727/[0 III]A5007; this combination is more sensitive 
to reddening, but has the benefit of tracing the ionization of a single element. Figure 4 shows 
that the filaments overlap the low-ionization extreme of LINER behavior, and in particular 
extend to very high [O l]/[0 III] ratios. We note further that in Figure 3 the measurements 
do not uniformly populate the LINER locus, but instead are clearly concentrated to rela- 
tively low values for each of the forbidden-line/recombination-line ratios. Another way of 
stating this result is that the nebulosity is not merely described by low ionization, as for 
LINERs, but is also characterized by low excitation. Since the forbidden fines arise from 
coUisional excitation, while the Balmer lines trace recombination, this result implicates a 
relatively low degree of heating per ionization in the cluster filaments. 

An interesting result illustrated by the 2-D spectra is that the distribution of points 
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shows a larger variation in the [O III]/H/3 ratio than in ratios employing the forbidden hnes 

of lower ionization, such as [N Il]/Ha, [S Il]/Ha, and [O I]/Hq;. The measurements include 
data for PA = 62° at the location of a giant H II region reported by Shields & Filippenko 
(1990a). The ratios for this locale are distinct from those elsewhere in the filaments, and 
contribute excursions to low values of [O I] /Ha and [S II] /Ha seen in Figures 36 and 3c. 
Removal of these points would increase the contrast between variation in [O III]/H/3 and in 
the other ratios. 



3.2. Coronal Emission and Cooling 

A set of weak lines that are of special interest are the coronal lines. These high-ionization 
coUisionally excited lines are emitted from gas at a temperature of ~ iq5.5-6.3 which is 
intermediate between that of the ICM, 10^~^ K, and the optical filaments, ~ 10^ K. The 
coronal lines trace gas in the process of cooling and their strength is a direct measure of 
the mass coohng rate, M (Hu, Cowie, & Wang 1985; Sarazin & Graney 1991). Detection of 
coronal lines may thus provide important diagnostics of the cooling-flow phenomenon. 

The wavelength range covered in our spectra includes several coronal lines, specif- 
ically [Ni XII]A4232, [Ar XIV]A4414, [Fe XIV]A5303, [Ca XV]A5445, [Fe X]A6374, and 
[Ni XII]A6702. None of these transitions was detected in our composite spectrum. Up- 
per limits to these lines were measured after subtracting a fit to the stellar continuum using 
a galaxy template of NGC 205, as described in §2. The location of [Fe X]A6374 coincides 
with the red wing of [O I]A6364 and consequently warranted special treatment. The fol- 
lowing procedure was used to remove the oxygen line. A stellar continuum was fitted over 
the wavelength range of 6250 to 6350 A as described in §2, simultaneously with a Gaussian 
to [O I]A6300. To remove [O I]A6364, we divided the template of [O I]A6300 by the line 
ratio of 3.0 dictated by atomic parameters (Osterbrock 1989), and subtracted the resulting 
profile and continuum fit from the composite spectrum. An upper limit for [Fe X]A6374 was 
derived from the residual. The upper limits on the coronal lines were corrected for extinction 
corresponding to Ay — 1.2 mag, the same as that applied to the detected hnes. The results 
are listed in Table 2. 

Utilizing theoretical models and our upper limits, we can calculate upper limits on the 
coohng rate of the hot ICM as measured in our apertures. The coronal fine expected to 
have the highest emissivity is [Fe X] A6374, for which Sarazin & Graney (1991) predict 
1.39 X 10^^ erg s~^ for every 1 M0 yr"-*^. For the one-pixel equivalent aperture corresponding 
to the composite spectrum, the upper limit to the reddening-corrected surface brightness of 
this line that we measure is ~ 2.7 x 10"^'' erg s~^ cm~^ arcsec"^. The composite spectrum 
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is averaged over regions that extend out to a radius of ~ 28". Assuming that the coronal 
emission is distribiitcd uniformly within this radius, our measurement suggests an upper 
limit for the luminosity of [Fe X] of ~ 4.0 x 10'^'' erg s~^, corresponding to a mass cooling 
rate of less than 2900 Mq yr~^. This value of M is well above the accretion rate derived 
from X-ray observations; upper limits to the other coronal features would lead to even less 
restrictive limits on M. These results underscore the practical difficulty in using the coronal 
fines to constrain M, even with the investment of substantial amounts of observing time. 

The strength of coronal cmisson relative to Ha is of interest if the nebular filaments are 
energetically coupled to the cooling ICM. A cold filament irradiated by an isobarically cooling 
gas reprocesses the ionizing radiation into Ha with an efficiency ~3% (Voit & Donahue 1990; 
Donahue & Voit 1991). The Ha strength would thus be 0.03 (5/2) {kT/i^mp)M, where nip 
is the proton mass, /i is the mean mass per particle in units of nip, T is the temperature of 
the hot gas, and k is the Boltzmann constant. In this picture, both the Ha and coronal-line 
luminosities scale with M, leading to well-defined predictions for their ratios, which we list 
in Table 3 (Voit, Donahue, & Slavin 1994; Sarazin & Graney 1991). 

The observed upper limits to the coronal-line strengths relative to Ha are larger than 
those of the predictions, thus leaving open the possibility that the filaments are ionized 
by the cooling ICM, a scenario we examine in more detail in Section 4. The same data 
set as employed here was used by Shields & Filippcnko (1992) to obtain a more restrictive 
limit on [Fe X]A6374, after binning the data on the basis of Ha surface brightness. The 
use of an intermediate surface brightness bin produced a somewhat higher S/N ratio in the 
continuum. Their upper limit to [Fe X]A6374/Ha was 1.4 x 10~^, adjusted to use the same 
upper-limit definition and reddening correction employed here. This value is ~ 36% less than 
the theoretical prediction, which may implicate other ionization sources for the filaments. 
However, this conclusion assumes that the cooling plasma irradiating the nebular gas lies 
within our measurement aperture, rather than on larger scales, and that the coronal feature 
displays the same velocity width as Ha. Similar observational upper limits for coronal lines 
from cooling-flow fllaments in other clusters have been reported by Donahue & Stocke (1994) 
and Yan & Cohen (1995). 

An additional complication to the use of the coronal lines as diagnostics of cooling 
stems from persistent uncertainties in the relevant atomic data (e.g., Mohan, Hibbcrt, & 
Kingston 1994; Pelan & Berrington 1995; Storey, Mason, & Saraph 1996; Ferguson, Korista, 
& Ferland 1997; Oliva 1997). As summarized by Voit et al. (1994) for [Fe XIV]A5303, 
alternative choices of atomic parameter values can lead to order-of-magnitudc differences 
in the predicted emissivities for a cooling ICM. As noted in Table 3, our upper limit to 
[Fe XIV]A5303/Ha is less than the value predicted by Voit et al. (1994) but larger than 
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that of Sarazin & Graney (1991). An unambiguous interpretation of coronal-line limits or 
detections will require further improvement in the atomic data. 

3.3. Ionization Mechanisms 

A question central to the understanding of cooling-flow filaments is the issue of what 
ionizes and powers these systems. A comparison of the recombination expected for mass 
deposition rates obtained from X-ray observations with the observed Ha line luminosities 
makes it clear that simple cooling of the gas with a transition through a nebular phase 
is insufficient to explain the energy requirements of the filaments in NGC 1275 and other 
objects (e.g., Heckman et al. 1989). A variety of mechanisms have been proposed for the 
ionization of cooling-fiow filaments. Below we compare the predictions of these scenarios 
with the results of our observational study, taking into account various physical constraints. 

3.3.1. Ionization by an AGN 

NGC 1275 shows evidence for a Seyfert-like nucleus (Seyfert 1943), as is clear from the 
emission-line spectrum (e.g., Filippenko & Sargent 1985) and nonstellar continuum from the 
central regions. If the extended nebulosity is photoionized by this continuum, the ionization 
state of the gas will reflect the ionization parameter, which we deflne here as the ratio of 
ionizing photon and electron densities at the face of the irradiated cloud. The fine ratio 
[S II]A6717/[S II]A6731 is sensitive to the electron density, but remains consistent with the 
low-density hmit {ue ^ 10^ cm~^; Osterbrock 1989) throughout the filaments. While the 
density may vary below this threshold within the measured regions, pressure arguments 
suggest that large fiuctuations are unlikely (e.g., Johnstone & Fabian 1988; Heckman et al. 
1989). If the nebular filaments are ionized by the AGN, we would then expect that the 
ionization parameter would largely vary in the same way as the ionizing fiux, proportional 
to r"^. The line ratio that is the best diagnostic of the ionization state is [O III]A5007/H/9, 
since this ratio involves a high- ionization line (e.g., Ferland & Netzer 1983; Ho et al. 1993). 
The filaments range from ~ 10" to ~ 30" in projected distance from the nucleus. If these 
structures lie in the plane of the sky, then the change in the ionization parameter due to 
geometric dilution of fiux would be approximately one order of magnitude, translating into 
a significant variation in [O IIl]/H/9. Values of this line ratio as a function of radius for each 
of the filaments arc shown in Figure 5. While negative gradients are seen in some places, 
the variations at large radial distance tend to be small and/or non-monotonic, a result that 
is in agreement with previous observations having more limited PA coverage (Johnstone & 
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Fabian 1988). 

The observed behavior may still be consistent with photoionization by a central source 
if the physical distance from the nucleus is not indicated by the projected distance. John- 
stone & Fabian (1988) constructed a detailed model to study the effects of geometry in the 
filaments of NGC 1275, and were able to reproduce the observed hne ratios as a function of 
projected distance, but suggested some disagreement with the observed surface brightness 
of the filaments. The present data show reasonable consistency between measured surface 
brightnesses and predictions for ionization by the nucleus. The flux in the ultraviolet con- 
tinuum measured by the International Ultraviolet Explorer (Kinney et al. 1991) and X-ray 
measurements from Einstein (Branduardi-Raymont et al. 1981) imply a 2-point spectral 
index between 2500 A (after dereddening) and 2 keV of aox = —1-4 (assuming fi, cx u"). Ex- 
trapolating the UV flux assuming a power law with this slope implies an ionizing luminosity 
of 2 X 10^^ photons s~^. Clouds with a density of rig = 100 cm~^ at a characteristic distance 
of ~ 7 kpc will be described by an ionization parameter of C/ ~ 10~^, consistent with the 
observed line ratios. These numbers also predict for ionization-bounded clouds an emission 
measure (EM) of ~ 500 cm^^ pc, in reasonable agreement (given the uncertainties in this 
calculation) with the observed (reddening-corrected) surface brightness of Ha of 2 x 10~^^ 
ergs s~^ cm~^ arcsec"^, corresponding to EM ?a 900 cm"^ pc. However, despite this nominal 
agreement, ionization by the AGN appears unlikely based on other line diagnostics that are 
independent of projection. 

A line of particular interest for this purpose is He II A4686 since it is a gauge of the 
hardness of the ionizing continuum. Figure 6 shows the composite filament spectrum for 
NGC 1275 in the region of He II A4686 together with the template continuum. The observed 
3(7 upper hmit for He II A4686/H/3 is 0.024, after dereddening. Photoionization simulations 
show that this value is sensitive to the ionization parameter for a given spectral shape. 
Based on the continuum and nebular parameters described above, we examined this issue 
using photoionization calculations obtained with CLOUDY, version 90.04 (Ferland 1996). 
The result is a prediction that He II A4686/H/3 ~ 0.10, which is four times our upper limit. 
Therefore, the nebular filaments in NGC 1275 are highly unlikely to be powered by the AGN. 

In addition to being a source of ionizing radiation, the AGN in NGC 1275 is a radio 
source (3C 84) with extended radio structures (Pedlar 1981 et al.; Pedlar et al. 1990), 
that may in principle power the optical nebulosity. The relation between the radio plasma 
and the nebular gas was discussed briefly by McNamara, O'Connell, & Sarazin (1996) in 
a multiwavelength study dealing with NGC 1275 and its environs. The extended radio 
emission roughly coincides with two cavities in the X-ray gas located on the northeast and 
southwest outer portions of the galaxy (Pedlar et al. 1990; McNamara et al. 1996). The 
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nebular emission seems to avoid the radio emission, which may be evidence for the radio 
lobes displacing the nebular gas. Another interpretation would be that the nebular gas, 
which has a higher density, may be confining the radio-emitting gas and channeling it to 
regions of lower density (McNamara et al. 1996). In any case, if an interaction with the 
radio lobes were responsible for powering the nebular emission, we might expect the optical 
filaments to be correlated spatially with the radio structures. This correspondence is not 
seen, but might indicate that the shocked gas remains in a hot phase. Recent observations 
with Chandra, however, show little evidence for X-ray emission by a hot shocked component 
in the regions in question (M. Wise, private communication); a similar situation has been 
reported for Hydra A (McNamara et al. 2000). As will be argued below, these results and 
additional spectroscopic constraints suggest that shock heating is unlikely in the filaments 
of NGC 1275. 



3.3.2. Shocks 

The above arguments point toward a distributed ionizing source as the underlying power 
source for the filaments in NGC 1275. This ionization mechanism could in principle be 

shocks, stars, and/or the cooling ICM itself. Several observational properties can be inter- 
preted as possible evidence of shocks in the filaments. The velocity field of the nebular gas 
shows variations of 100 — 300 km s~^ (e.g., Burbidge & Burbidge 1965; Heckman et al. 1989). 
Moreover, the line widths are greater than the thermal width. Line emission from gas at a 
temperature of 10,000 K will exhibit a full width at half maximum (FWHM) of ~ 20 km s^^ 
for hydrogen, which can be compared with observed FWHM values of 100 — 300 km in 
the spatially-resolved spectra, after removal of the instrumental width. 

Several possible mechanisms are available for driving shocks, including interaction of 
the radio plasma from the AGN with surrounding gas (§3.3.1), pressure imbalances resulting 
from thermally unstable coohng in the cluster center (Cowie, Fabian, & Nulsen 1980), and 
accretion infall onto the central galaxy. For shocks driven by radio ejecta, Bicknell, Dopita, 
& O'Dea (1997) have shown that the luminosity in Ha is of order 0.1% of the energy flux into 
the radio plasma, with the exact value dependent on the gas density distribution and shock 
velocity. For NGC 1275, Pedlar et al. (1990) estimate that ~ 10'^^ erg is injected into 
the radio halo; this value is inversely related to the age of the halo, which Pedlar et al. take 
to be 10^ years or less. The total Ha luminosity of the NGC 1275 filaments is ~ 2.1 x 10^^ 
erg s~^ (Heckman et al. 1989; Caulet et al. 1992). It appears doubtful that the radio source 
is sufficient to power the observed line emission, unless the radio halo is exceedingly young 
or the energy transfer to optical emission is unusually efficient. 
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Efficient cooling of matter near the cluster center can lead to regions that are under- 
pressured relative to their surroundings, and hence will be subjected to repressurizing shocks 
(Cowie, Fabian, & Nulsen 1980). The power delivered by such shocks results from the ther- 
mal and mechanical energy of the ICM, so that the maximum power that be can extracted 
from the cooling ICM is (5/2) {kT/fimp)M ^ 10^^ erg s~^. Only about 3% of this energy 
is processed into Ha, requiring that the cooling clouds be shocked multiple times to repro- 
duce the observed emission (the "i?rec problem"; see Heckman et al. 1989 for a discussion). 
The apparent survival of grains in the filaments of other cooling-fiow galaxies suggests that 
passage of repeated shocks through the nebular gas is unlikely (Donahue & Voit 1993). 

If the bulk kinetic energy of infalling matter dominates the energy inflow to the central 
region, the power available for driving shocks is of order Mv"^, where Vi is the characteristic 
infall speed. If Vi is as high as the cluster velocity dispersion (~ 1000 km s~^; Kent & Sargent 
1983), the resulting rate of energy transfer is ~ 10^^ erg s""^, which is sufficient to power the 
filament emission. However, if Vi is closer to the filament velocity dispersion, the available 
power is an order of magnitude less, and probably insufficient to generate the observed Ha 
luminosity. 

One general characteristic of shocks is that line-intensity ratios of emitted radiation 
refiect the kinetic energy flux entering the shock front, and hence are a function of shock 
velocity (e.g., ShuU & McKee 1979; Dopita & Sutherland 1995, Binette, Dopita, & Tuohy 
1985). If the observed line width reflects the characteristic shock velocity, then we would 
expect the line ratios to correlate with line width; such behavior is in fact seen in large-scale 
outflows from starburst galaxies (e.g., Heckman et al. 1999). A plot of log([N II]A6583/Ha) 
versus the FWHM of Ha from the four PAs is shown in Figure 7, together with model 
predictions, which exhibit a positive trend in most cases. The range of observed FWHM 
is 100 — 300 km s~^, while [N Il]/Ha is clustered around 0.9. A shallow but statistically 
significant correlation (Spearman non-parametric significance of ~ 5 x 10~^) is present, with 
larger line ratios associated with higher velocity width. This trend is due to combining values 
from all the PAs; PA = 283° and PA = 313° have slightly higher line ratios and FWHM 
than PA = 62° and PA = 77°. Significant correlations are not seen in either the individual 
PAs or in PA62/PA77 and PA283/PA313 groupings, making shocks as heating sources in the 
filaments unlikely. Similar behavior is seen in the [O III]/H/3 ratios, although no correlation 
or PA-dependent behavior is evident in [S II] /Ha or [O I] /Ha. 

A further piece of evidence which argues against shocks comes from the upper limit 
on electron temperature as determined from the [O III] lines (Table 1). Figure 8 shows 
the continuum around the expected location of [O III]A4363. An upper limit on the tem- 
perature as measured from [O IIl]A4363/[0 III]AA4959, 5007 <0.047 restricts the electron 
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temperature to be less than 33,000 K (Osterbrock 1989; Shaw & Dufour 1995). Shock mod- 
els predict [O IIl]A4363/[0 III]AA4959, 5007 to be between 0.01 and 0.07 (ShuU & McKee 
1979; Raymond 1979; Binette et al. 1985; Dopita & Sutherland 1995), corresponding to 
electron temperatures between 13,000 K and 54,000 K. Shock models with velocities less 
than ~ 300 km s^^ predict an [O III] line ratio consistent with our upper limit, but have 
problems in simultaneously matching the high- and low-ionization forbidden-line strengths. 
These models also consistently underpredict He I A5876, except at high velocities, where the 
predicted He II A4686 exceeds our upper limit (Raymond 1979; Binette et al. 1985; Dopita & 
Sutherland 1995). Therefore, we conclude that shock heating in the filaments of NGC 1275 
is a minor contributor to the energetics of the nebular gas. 



3.3.3. Stellar Photoionization 

Cooling-flow galaxies commonly exhibit an optical continuum that is bluer than typical 
elliptical galaxies (e.g., Johnstone, Fabian, & Nulsen 1987; McNamara & O'Connell 1989; 
Allen 1995; Cardiel, Gorgas, & Aragon-Salamanca 1995, 1998). Some degree of correlation 
exists between such blue continua and nebular emission, prompting suggestions that both 
may trace the influence of an underlying young stellar population. While nebular emission 
(bremsstrahlung, free-bound, and 2-photon) may contribute to such a continuum, a quan- 
titative calculation indicates that the continuum emission accompanying the observed H/? 
luminosity would be 2 orders of magnitude too weak to account for the observed light. Stars 
are thus the likely source of the blue continua. As possible ionizing agents for the nebular gas, 
stars also have the added appeal of producing negligible He II nebular emission, consistent 
with the observations of NGC 1275 and other cooling-flow galaxies. 

Early studies of the large-scale nebular morphology of NGC 1275 interpreted this object 
and its system of radial filaments as the result of a large-scale outfiow; indeed, it was called 
an "exploding galaxy" (Burbidge & Burbidge 1965). Modern studies of another object given 
this label, M82, demonstrate that it is in a sense "exploding" because of its large-scale 
outfiow driven by star formation (e.g., McCarthy et al. 1987; Shopbell & Bland-Hawthorn 
1998). However, such vigorous star formation is not apparent in NGC 1275. 

NGC 1275 does show weaker indications of star formation within the last ~ 10^ — 10^ 
yrs on a variety of scales (Wirth et al. 1983; Holtzman et al. 1992; McNamara et al. 1996). 
Optical measures of the young population have been used to estimate roughly the associated 
ionizing luminosity, yielding a result that is consistent with the requirements of powering 
the filaments (Wirth et al. 1983; Romanishin 1987). However, there is not a strong spatial 
correspondence between the regions of recent star formation and the nebular gas. As noted 



-14- 



previously, the nebular spectrum also differs markedly from typical H II region spectra; in 
the one instance where recent star formation is clearly evident within a filament (at PA = 
62°), the nebular line ratios are distinct from the remainder of the filament system (Shields 
& Filippenko 1990a). Photoionization models suggest that hot stars can generate LINER- 
like emission (Filippenko & Terlevich 1992; Shields 1992), but these scenarios generally have 
difficulty producing [O I] and [N I] emission as strong as seen here, unless the gas contains 
a very high density component (~ 10^ cm~^). 

Stars may nonetheless play a role in powering the filaments if they generate a radiation 
field that ionizes gas on scales exceeding that of localized H II regions. A similar scenario has 
recently gained favor for explaining the ionization of the diffuse ionized gas (DIG) commonly 
seen in spiral galaxies (e.g., Howk & Sembach 1999; Rand 1998; Martin 1997). However, the 
forbidden-hne strengths seen in the DIG are characteristically weaker than in the NGC 1275 
filaments; typical DIG line ratios are [O I]A6300/Hq; ^ 0.03, [N II]A6583/Hq; ^ 0.5, He I 
A5876/Ha ^ 0.02, and [O III]A5007/Ha ^ 0.1 (Galarza, Waherbos, & Braun 1999; Sembach 
et al. 2000), which contrasts markedly with the values in Table 1. Another aspect which 
makes the cooling-fiow filaments of NGC 1275 distinct from the DIG is the emission measure. 
For the DIG, EM ^ 1 - 10 cm^^ pc (Reynolds 1985; Martin & Kennicutt 1997; Hoopes, 
Walterbos, & Rand 1999), and can reach values as high as 100 cm^® pc (Galarza et al. 1999; 
Hoopes et al. 1999). In contrast, the average EM in NGC 1275 for the regions studied here 
is ~ 900 cm~^ pc. 

Given the arguments for and against stellar photoionization, it may be the case that 
this excitation process, if present, is at work in combination with another ionization source, 
as has been pointed out by Voit & Donahue (1997) for Abell 2597. The evidence from 
NGC 1275, however, does not form a strong case for stellar photoionization as the dominant 
heating mechanism in the filaments. 

3.3.4- Photoionization by the ICM 

The hot ICM surrounding NGC 1275 can be the ionization source of the filaments if 
its energy is channeled efficiently into the nebulosity (Voit & Donahue 1990; Begelman & 
Fabian 1990). In one such scenario, known as self- irradiation, the hot ICM ionizes nebulae 
that have already cooled and condensed from the ICM itself (Donahue & Voit 1991; Voit & 
Donahue 1990; Voit et al. 1994). Another way to channel ICM energy into the nebulae is 
the formation of mixing layers between cold clouds (Te ^ lO"^ K) embedded in a hot ambient 
medium {Tg ^ 10^ K) (Begelman & Fabian 1990; Crawford & Fabian 1992). Line emission 
comes from the mixing layer, which is described by an intermediate temperature (Te a; 10^'^ 
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K), and from the cool cloud that is now being irradiated by the mixing layer. Crawford & 
Fabian (1992) were able to reproduce the line ratios of Abell 496, a cluster that displays Type 
I filaments, with a mixing-layer prescription. However, in order to reproduce the line ratios 
of the Type 11 cluster Abell 2597, these authors found it necessary to invoke an additional 
energetic contribution from shocks. 

The observed line ratios in NGC 1275 are close to values predicted by one of the self- 
irradiation models presented by Voit et al. (1994). In that model (D39T67), the ICM is 
cooling from a temperature of 5x10^ K, the ionization parameter is 10^^'^, and the column 
density of irradiated H I is 10^^'^ cm~^. Their predicted He II A4686/H/9 ratio is almost a 
factor of 3 larger than our upper limit. This comparison shows that self-irradiation is not 
the ionization mechanism because it provides a continuum which is harder than that inferred 
from the observed line ratios. A cooling-flow system that shows similar filaments, in terms 
of line ratios, to NGC 1275 is Abell 2597 (Heckman et al. 1989; Voit & Donahue 1997). Voit 
& Donahue (1997) similarly ruled out a hard ionizing continuum as the energy source in the 
filaments of Abell 2597 on the basis of an upper limit of 0.02 on He II A4686/H/?, which is 
again a factor of ~ 3 below expectations for self-irradiation. 



3.3.5. Composite Ionization Models 

None of the ionization mechanisms considered individually in the preceding sections is 
entirely satisfactory in accounting for the emission properties and energetics of the cooling- 
fiow filaments in NGC 1275. Given the compelling evidence for the presence of young 
stars, it is reasonable to explore models that are composites of stellar photoionization and 
some secondary heating mechanism. The study by Voit & Donahue (1997) of Abell 2597 
is illuminating in this regard. Voit & Donahue (1997) invoked stellar photoionization as 
the ionizing mechanism for the optical filaments of Abell 2597 after placing strong upper 
hmits on [O III]A4363 and He II A4686, a result that ruled out shock heating and self- 
irradiation by the ICM. However, they found that the electron temperature predicted for 
stellar photoionization was lower than what they inferred from their observations of the 
[O II] and [S II] lines. For this reason, they suggested that a heating source in addition to 
hot stars may be important. The secondary source they considered was the ICM; it can 
heat the filaments via electron conduction (e.g.. Sparks 1992), acoustic-wave heating (e.g., 
Pringle 1989), and MHD-wavc heating (e.g., Friaca et al. 1997). Voit & Donahue (1997) 
argued that the resulting energy fiux incident on the nebulae can become greater than that 
of photoelectric heating if the nebulae are characterized by low ionization, a case which 
holds in cooling-fiow filaments. Unfortunately, detailed theoretical predictions appropriate 
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for comparison with observations have not yet been generated for these scenarios. A related 
case where the heating rates due to the dissipation of tTirbulence have been calcTilatcd and 
found to be comparable to the radiative cooling rates is in the DIG of the Milky Way (Minter 
& Spangler 1997). However, detailed observational parameters (such as the amplitude of 
density fluctuations, and turbulence length scales) needed to calculate the heating rate from 
this mechanism arc unavailable for cooling-flow fllaments. 

We investigated an alternative scenario in which the observed emission arises from a 
combination of clouds subject to stellar photoionization and irradiation by the ICM. In this 
simplified model, we assumed that within a given measurement aperture, we sec a mix of 
clouds ionized by stars and by radiation from the ICM. In this picture we observe a weighted 
average of the spectrum of the two types of cloud. A fraction of the observed fine emission 
comes from regions ionized by stars, while another fraction comes from nebulae irradiated 
by the ICM. This weighted average translates into the corresponding fractions in terms of 
line ratios (with all ratios normalized to H/3). We used line ratios from the self- irradiation 
models of Voit et al. (1994), and representative values for intermediate- mctallicity H II 
regions (Veilleux & Osterbrock 1987). Once the relative contributions of stars and the ICM 
were determined from He II A4686/H/9, attention was turned to reproducing [O I]A6300/Ha 
as well as the other line ratios. [O I]A6300/Hq; is a good diagnostic because it is very weak in 
H II regions. Including more stars pushed [O I]A6300/Hq! below observed values, while more 
self- irradiation would push He II A4686/H/3 above what is observed. Ultimately we were not 
able to find a composite model that fits all the line ratios simultaneously. Combinations 
satisfying the helium requirement, with marginal agreement for the other line ratios, failed 
to match the oxygen constraint, and vice versa. Composite models can be found that satisfy 
both [O I] /Ha and He II/H/?, but these scenarios proved inconsistent with the other line 
ratios. The reason for this failure to come up with a successful hybrid model may be due to 
our simplistic approach, which neglects clouds that sec ionizing radiation from both sources 
simultaneously, or more likely suggests that other ionization mechanisms are at work. 



3.3.6. Attenuated Cluster Radiation 

If the filaments in NGC 1275 and other cooling-fiow galaxies are photoionized, hints 
about the nature of the ionizing continuum come from the strengths of the emission lines. 
Relatively strong forbidden lines, the absence of He II A4686, and the presence of He I A5876 
suggest that the ionizing continuum cuts off at energies somewhere between 24 and 54 eV. 
An ionizing source generating a continuum with a natural cutoff at 54 eV is normal 0-type 
stars. For this reason we tried in the previous subsection to construct a composite model of 
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stellar photoionization and self-irradiation by the ICM, but our attempts were unsuccessful 
in terms of matching the observed line ratios. 

The cooling ICM can produce a range of ionizing continua, depending on its initial 
temperature. Self-irradiation models (Voit et al. 1994) that predict a He II A4686 value 
similar to what is observed fail at reproducing the other line ratios, while self-irradiation 
models that agree with most of the observed line ratios have high He II A4686 /H/9. One way 
to construct an appropriate ionizing spectrum with the shape dictated by the observations 
is to introduce a cutoff at the helium edge of 54 cV by attenuating the spectrum by an 
ionized screen. The filtering medium will imprint a sharp dip at 54 cV while leaving the flux 
of photons with energies between 24 and 54 eV nearly unchanged. Observational evidence 
for such a screen may come from X-ray observations of Abell 426 that suggest the presence 
of absorbing matter with a column density of Nh ~ 10^^ cm~^, in excess of the Galactic 
foreground component; the ionization state of this medium is not well understood (e.g., 
Allen, Fabian, & Edge 1992; Arnaud & Mushotzky 1998). 

We therefore constructed models in which the ionizing continuum from the cooling ICM 
is attenuated through a screen before irradiating the filaments. This scenario differs from the 
traditional self-irradiation models of Voit & Donahue (1990) in the sense that the filaments 
and the screen need not be the by-products of the cooling flow. The models discussed here 
also differ from the mixing layers picture of Crawford & Fabian (1992) in terms of the origin 
of the ionizing continuum. In the mixing layers model, the ionizing spectrum comes from gas 
at a temperature intermediate between the temperature of the ICM and that of the optical 
filaments. 

To test the photoionization plus screen scenario, we carried out photoionization calcu- 
lations using CLOUDY, version 90.04 (Ferland 1996). We assumed plane-parallel geometry 
throughout the calculations, and a mctallicity of 0.4 times solar values, consistent with X-ray 
studies (e.g., Allen et al. 1992). The fiux originating from the ICM was scaled to match the 
soft X-ray central surface brightness reported by Mohr, Mathiesen, & Evrard (1999), and 
propagated consistently through the calculation. We found good results assuming a tem- 
perature for the dominant ionizing plasma of 10^ K. While lower than that of the majority 
of the ICM, this value is in general accord with the high degree of cooling occurring in the 
Perseus cluster core, and with spectroscopic evidence for cooler phases (e.g., Mushotzky & 
Szymkowiak 1988) in this region. 

In order to have the appropriate screening properties, the ionization state of the filtering 
medium should be such that the optical depth at 54 eV is large enough to suppress He II 
A4686 in the filaments, while the total hydrogen column density stays less than, or com- 
parable to, the column density inferred from X-ray observations. Moreover, the ionization 
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parameter of the screen should be high enough to produce minimal optical forbidden-line 
emission, since the screen is evidently not seen in sTich lines. Therefore, we chose the hydro- 
gen density in the screen to be 1 cm~^, implying an ionization parameter of ~ 10^^'^, and 
stopped the calculation when the optical depth at 54 eV became 10. The resulting hydrogen 
column density in the screen is ~6xl0^^ cm~^. For the filament, the forbidden-line spec- 
trum suggests an ionization parameter U ~ 10~^'^ (§3.3.1). This value emerges naturally 
from the propagated radiation field and a choice of density of ~ 100 cm~^, consistent with 
the observed [S II]A6717/A6731 line ratio. 

Resulting predictions for the emission-line spectrum emergent from a filament are shown 
in Table 3. All the line ratios were satisfactorily reproduced, with the exception of the [S II] 
lines, which are too strong by a factor of ~ 1.6. Figure 9 shows the continuum incident on 
the screen, and the continuum exiting the screen and incident on the filament. Notice the 
strong ionized helium edge at 54 eV. Since the ionized screen is photoionized, it will also 
emit recombination radiation, especially in He II A4686 (He II Paa). In a sense, then, this 
model removes the difficulty of excess He II emission in the filaments, but probably creates 
a new observational inconsistency in that the screen should be a strong source of observable 
recombination emission, which is not apparently seen. 

He II A4686 emission would be reduced by a factor of ~ 2 if recombination in the screen 
were described by Case A conditions rather than Case B. However, the line center optical 
depth for He II Ly7, which upon reabsorption can branch to emit the Pa« line, remains ^ 1 
for the column density and temperature describing our model screen; additional turbulence 
of ~ lO'' — 10^ km s~^ within the screen would be required to yield Case A conditions. Case 
A recombination thus does not appear to be a plausible mechanism for reducing emission in 
He II A4686. 



3.3.7. Diffuse Extreme Ultraviolet Radiation 

Observations with the Extreme Ultraviolet Explorer (EUVE) indicate the presence of 
excess soft X-ray/EUV emission in clusters of galaxies (Lieu et al. 1996; Lieu, Bonamente, 
& Mittaz 1999). The existence and properties of this emission are still highly debated (e.g., 
Arabadjis & Bregman 1999). If this emission component is real, an interesting question is 
whether it plays a significant role in the energetics of the nebular filaments. 

A difficulty with using the EUV observations directly stems from the fact that EUVE 
data (Lieu et al. 1996) cover only a limited bandpass, making any attempts at determining 
the total ionizing luminosity dependent upon the shape of the assumed continuum, a property 
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that has to be determined by a model for the origin of the EUV excess. According to one 
interpretation, inverse- Compton scattering of the cosmic microwave background radiation 
off cosmic-ray electrons is responsible for the EUV continuum, which is then described by 
a two-point spectral index between 70 and 912 A of asuv = —2/3 (Sarazin & Lieu 1998; 
Sarazin 1999). With this spectral energy distribution, the ratio of photon surface brightness 
between 70 A and 912 A to that between 70 A and 160 A, the approximate range covered 
by the EUVE DS Lexan detector, is about 6.2. 

No EUVE observations have been made of NGC 1275; the large Galactic H I column 
density toward this source (1.4 x 10^^ cm~^; Lockman & Savage 1995) makes its detection dif- 
ficult or impossible at EUV wavelengths. In order to obtain an order-of-magnitude estimate 
of the relevant photon flux in the EUVE bandpass, we used observations of Abell 1795 (Lieu 
et al. 1999). Abell 1795 is a cooling-flow cluster whose central galaxy displays bright optical 
filaments which, hke those in NGC 1275, were classified spectroscopically by Heckman et al. 
(1989) as Type II; the M values inferred from X-ray observations agree to within a factor 
of ~ 2 for the two clusters. Assuming an effective area of ~ 15 cm^ for EUVE with the 
DS Lexan detector (Miller 1997), the central surface brightness for Abell 1795 is 2.0x10^^ 
photons s^^ cm^^ arcmin"^ (Lieu et al. 1999). Scaling this value to obtain the ionizing (70 A 
- 912 A) surface brightness results in a value of ~ 1.2 x 10~^ photons s~^ cm~^ arcmin"^. 
If we assume that the filaments can be treated as ionization-bounded slabs located at the 
center of the EUV-emitting medium, the ionizing photon fiux striking the optical gas is 
~ 4.5 X 10^ photons cm~^, implying EM !v 0.06 cm~^ pc in the photoionized slab. This 
value differs substantially from the emission measure calculated above from our Ha surface 
brightness for NGC 1275, ~ 900 cm~^ pc, which is also comparable to the results obtained 
for other systems (Voit & Donahue 1997). We conclude that the EUV emission, if present 
in Abell 426, is a highly unlikely ionizing mechanism for the optical filaments in NGC 1275. 

3.3.8. Magnetic Reconnection 

Another source of energy in galaxy clusters is the intergalactic magnetic field. Magnetic 
reconnection of field lines in clusters can release energy in a manner similar to that of solar 
flares (e.g., Heckman et al. 1989; Soker & Sarazin 1990). For magnetic reconnection to be 
effective in releasing energy, the magnetic field must be strong and tangled. Reconnection 
of magnetic flux loops will liberate energy which has to be included in the energy budget of 
the cooling-flow fllaments (e.g. Soker & Sarazin 1990; Jafelice & Friaca 1996). 

Cooling flows present in galaxy clusters may amplify the magnetic fleld energy density 
to equipartition values with the gas pressure within the inner tens of kpc due to compression 
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of the field lines that are frozen-in to the plasma (e.g., Soker & Sarazin 1990). It is interesting 
to note that the observed optical filaments in cooling-flow galaxies occur over the same length 
scale predicted for equipartition between magnetic and gas pressures, for theoretical cooling- 
fiow models. The resulting strong magnetic fields are expected to produce large Faraday 
rotations in the cores of cooling-flow clusters, in agreement with radio measurements for 
many sources (Ge & Owen 1993, 1994). As a result of this amplification, one would expect 
clusters with larger coohng rates to have larger Faraday rotation measures. 

The correlation between cooling rates and Faraday rotation measure seems to be borne 
out by observations. Taylor, Barton, & Gc (1994) compiled a list of radio sources at the 
centers of galaxy clusters, and showed that such a trend exists. Notable exceptions, however, 
are NGC 1275 and PKS 0745-191, which have strong coohng flows and extended optical 
nebulae, but no measurable polarization (Ge & Roberts 1993; Taylor et al. 1994). One 
possible interpretation of this finding is that a magnetized screen with complex structure is 
strongly depolarizing the radio signal (Ge & Roberts 1993; Taylor et al. 1994). The screen 
may be the optical filaments themselves. Several pieces of observational evidence indicate 
that this might be the case. The Ha and radio emission are cospatial in projection in 
PKS 0745-191 and in the core of NGC 1275 (Heckman et al. 1989; Ge & Roberts 1993). Van 
Breugel, Heckman, & Miley (1984) noted that there is a spatial inverse correlation between 
radio polarization and Ha surface brightness in Abell 1795, a result later confirmed by Ge 
& Owen (1993) at higher spatial resolution. Obtaining more radio polarimetry data with 
high spatial resolution for cooling-flow clusters will be very helpful in testing this inverse 
correlation. A detailed rotation measure map will make it possible to look for a relation 
between the Ha luminosity of a filament and the degree of depolarization it imprints on the 
radio emission cospatial with it. 

The above discussion provides circumstantial evidence for the importance of magnetic 
fields in understanding cooling-flow fllaments. If Faraday depolarization in the optical flla- 
ments can be taken as evidence for the fllaments being magnetized with a complex geometry, 
then magnetic field lines may be reconnecting in the nebular gas, and one would expect strong 
depolarization to be associated with strong optical emission. In terms of energetics, Jafelice 
& Friaca (1996) argued, on the basis of theoretical models, that magnetic reconnection is 
capable of powering the lower luminosity Type I filaments, but falls short of sustaining the 
Type II systems, such as NGC 1275, AbeU 1795, and PKS 0745-191. Additional high spatial 
resolution radio polarimetry for cooling-fiow galaxies is desirable in order to carry out a 
comparative study of the magnetic properties of the two source types. A further interesting 
experiment would be to obtain sTich data for the central galaxy in Sersic 159-03 (S 1101), 
which displays filaments with emission-line ratios encompassing both Type I and Type II 
values (Crawford & Fabian 1992). 
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4. Conclusions 

In this paper we have presented the results of a detailed optical spectroscopic study of 
the optical filaments around NGC 1275, the central galaxy in the Perseus cluster of galaxies. 
The database was extensive enough to yield tight upper limits on weak lines of interest for 
constraining the ionization and thermal state of the nebular gas. Line-ratio diagrams reveal 
that the nebulosity has an unusual combination of emission-line strengths reflecting both 
low ionization and low excitation. 

The spectroscopic results were used to gauge the viability of several scenarios proposed 
for powering cooling-fiow nebulosity. Mechanisms examined include shocks, ionization by an 
AGN central source or extended radio structures, stellar photoionization, and energy transfer 
by various means from the ICM. None of these models is entirely satisfactory in accounting 
for the energetic and/or emission-line properties of the filaments. Radio studies of Faraday 
rotation for cooling-fiow galaxies provide tantalizing suggestions that magnetic reconnection 
plays a significant role in the filament energetics. Additional polarization studies of these 
sources may provide a clearer picture as to the importance of magnetic fields in these systems. 

Support for this work was provided in part by a grant from the California Space Institute 
to J.C.S., by NSF grants AST-8957063 and AST-9417213 to A.V.F., and by NASA grant 
NAG5-3556 to A.V.F. We thank the Mount Hamilton staff of Lick Observatory for their 
friendly and capable assistance over the multiple observing runs required for this project. 
At the time when the observations were obtained. Lick Observatory received partial funding 
from the NSF through Core Block Grant AST-8614510 to the University of California. We 
also thank Ray Bertram and Rick Pogge for assistance at the Perkins 1.8-m telescope. We 
gratefully acknowledge Gary Ferland for making CLOUDY available for this study, and Mike 
Wise for communicating results from Chandra observations in advance of publication. 

REFERENCES 

Arabadjis, J. S., & Bregman, J. N. 1999, ApJ, 514, 607 
Allen, S. W. 1995, MNRAS, 276, 947 

Allen, S. W., Fabian, A. C, Johnstone, R. M., Arnaud, K. E., & Nulsen, P. E. J. 2000, 
MNRAS, submitted ( |astro-ph/9910188|) 

Allen, S. W., Fabian, A. C, Nulsen, P. E. J., & Edge, A.C. 1992, MNRAS, 254, 51 

Arnaud, K. A., & Mushotzky, R. F. 1998, ApJ, 501, 119 



-22- 



Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Begelman, M. C, & Fabian, A. C. 1990, MNRAS, 244,26P 
Bicknell, G. V., Dopita, M. A., & O'Dea, C. P. O. 1997, ApJ, 485, 112 
Binettc, L., Dopita, M. A., & Tuohy, I. R. 1985, ApJ, 297, 476 

Branduardi-Raymont, G., Fabricant, D., Feigelson, E., Gorenstein, J., Grindlay, J., Soltan, 

A., & Zamorani, G. 1981, ApJ, 248, 55 

Burbidge, E. M., & Burbidge, G. R. 1965, ApJ, 142, 1351 

Burstein, D., & Heiles, C. 1984, ApJS, 54, 33 

Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 

Cardiel, N., Gorgas, J., & Aragon-Salamanca, A. 1995, MNRAS, 277, 502 

Cardiel, N., Gorgas, J., & Aragon-Salamanca, A. 1998, MNRAS, 298, 977 

Caulet, A., et al. 1992, ApJ, 388, 301 

Cowie, L. L., Fabian, A. C., & Nulsen, P. E. J. 1980, MNRAS, 191, 399 

Crawford, C, & Fabian, A. C. 1992, MNRAS, 259, 265 

Donahue, M., & Stocke, J. T. 1994, ApJ, 422, 459 

Donahue, M., & Voit, G. M. 1991, ApJ, 381, 361 

Donahue, M., & Voit, G. M. 1993, ApJ, 414, L17 

Dopita, M. A., & Sutherland, R. S. 1995, ApJS, 102, 161 

Fabian, A. C, Cowie, L. L., & Grindlay, J. 1981, ApJ, 248, 47 

Fabian, A. C, Nulsen, P. E. J., & Canizares, C. R. 1984, Nature, 310, 733 

Ferguson, J. W., Korista, K. T., & Ferland, G. 1997, ApJS, 110, 287 

Ferland, G. J. 1996, Hazy, a Brief Introduction to CLOUDY, University of Kentucky De- 
partment of Physics and Astronomy Internal Report 

Ferland, G. J., & Netzer, H. 1983, ApJ, 264, 105 

Filippcnko, A. V., & W. L. W. Sargent 1985, ApJS, 85, 503 

Filippenko, A. V., & Terlevich, R. 1992, ApJ, 397, L79 

Friaca, A. C. E., Gonvalves, D. R., Jafelice, L. C, Jatenco-Pereira, V., & Opher, R. 1997, 
MNRAS, 324, 449 

Galarza, V. C, Walterbos, R. A. M., & Braun, R. 1999, AJ, 118, 2775 

Ge, J. P., & Owen, F. N. 1993, AJ, 105, 778 



-23- 



Ge, J. P., & Roberts, J. D. 1993, BAAS, 183, 53.10 
Heckman, T. M. 1980, A&A, 87, 152 

Heckman, T. M., Armus, L., Weaver, K. A., & Wang, J. 1999, ApJ, 517, 130 

Heckman, T. M., Baum, S. A., van Breugel, W. J. M., & McCarthy, P. J. 1989, ApJ, 338, 
48 

Held, E. v., Mould, J. R., & dc Zccuw, P. T. 1990, AJ, 100, 415 
Holtzman, E. J., et al. 1992, AJ, 103, 691 

Hoopes, C. G., Walterbos, R. A. M., & Rand, R. 1999, ApJ, 522, 669 

Howk, J. C., & Sembach, K. R. 1999, ApJ, 523, L141 

Hu, M. E., Cowie, L. L., & Wang, Z. 1985, ApJS, 59, 447 

Hummer, D. G., & Storey P. J. 1987, MNRAS, 224, 801 

Jafelice, L. C, & Friace, A. C. S. 1996, MNRAS, 280, 438 

Johnstone, R. M., & Fabian, A. C. 1988, MNRAS, 233, 581 

Johnstone, R. M., Fabian, A. C., & Nulsen, P. E. J. 1987, MNRAS, 224, 75 

Kent, S. M., & Sargent, W. L. W. 1979, ApJ, 230, 667 

Kent, S. M., & Sargent, W. L. W. 1983, AJ, 88, 697 

Kinney, A. L., Bohlin, R. C., Blades, J. C., & York, D. G. 1991, ApJS, 75, 645 

Kriss, G. 1994, in ASP Conf. Ser. 61, Astronomical Data Analysis Software and Systems 
HI, ed. D. R. Crabtree, R. J. Hanisch, & J. Barnes (San Francisco: ASP), 437 

Lieu, R., Bonamente, M., & Mittaz, J. P. D. 1999, ApJ, 517, L91 

Lieu, R., Mittaz, J. P. D., Bowyer, S., Lockman, F., & Hwang, C. 1996, ApJ, 458, L5 

Lockman, F. J., & Savage, B. D. 1995, ApJS, 97, 1 

Martin, C. L. 1997, ApJ, 491, 561 

Martin, C. L., & Kennicutt, R. C. Jr. 1997, ApJ, 483, 698 
McNamara, B. R., & O'Connell, R. W. 1989, AJ, 98, 2018 
McNamara, B. R., O'Connell, R. W., & Sarazin, C. L. 1996, AJ, 112, 91 
McNamara, B. R., et al. 2000, ApJ, 534, L135 

Miller, A. 1997, EUVE Guest Observer Program Handbook, 4th Ed. (Berkeley, CA: EGO 
Center) 



-24- 



Miller, J. S., & Stone, R. P. S. 1987, Lick Observatory Technical Report No. 48, The CCD 
Cassegrain Spectrograph at the Shane Reflector (Santa Cruz: Lick Observatory) 

Minkowski, B. R. 1957, in lAU Symp. 4, Radio Astronomy, ed. H. C. van der Hulst 
(Cambridge: Cambridge Univ. Press), 107 

Minter, T., & Spangler, S. R. 1997, ApJ, 485, 182 

Mohan, M., Hibbert, A., & Kingston, A. E. 1994, ApJ, 434, 389 

Mohr, J. J., Mathiesen, B., & Evrard, A. E. 1999, ApJ, 517 

Mushotzky, R. P., & Szymokowiak, A. E. 1988, in Cooling Flows in Clusters and Galaxies, 
ed. A. C. Fabian, (Dordrecht: Kluwer), 53 

Nelson, C. H., & Whittle, M. 1995, ApJS, 99, 67 

Oliva, E. 1997, in ASP Conf. Scr. 113, Emission Lines in Active Galaxies: New Methods 
and Techniques, ed. B. M. Peterson, F. Cheng, & A. S. Wilson (San Francisco: ASP), 

228 

Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill 
Valley, CA: University Science Books) 

Pedlar, A., Booler, V., & Davies, R. V. 1983, MNRAS, 203, 667 

Pedlar, A., Ghantaure, H. S., Davies, R. V., Harrison, B. A., Perley, R., Crane, P. C, & 
Unger, S. W. 1990, MNRAS, 246, 477 

Pelan, J., & Berrington, K. A. 1995, A&A, 110, 209 

Pringle, J. E. 1989, MNRAS, 239, 479 

Rand, R. 1998, ApJ, 501, 137 

Raymond, J. C. 1979, ApJS, 39, 1 

Reynolds, R. J. 1985, ApJ, 294, 256 

Romanishin, W. 1987, ApJ, 323, L113 

Rubin, V. C, Ford, W. K., Peterson, C. J., & Lynds, C. R. 1978, ApJS, 37, 235 

Sarazin, C. L. 1999, ApJ, 520, 529 

Sarazin, C. L., & Lieu, R. 1998, ApJ, 494, L177 

Sarazin, C. L., & Graney, C. M. 1991, ApJ, 375, 532 

Scmbach, K. R., Howk, J. C., Ryans, R. S. I., & Keenan, F. P. 2000, ApJ, 528, 310 

Seyfert, C. A., 1943, ApJ, 97, 28 

Shaw, R. A., & Dufour, R. J. 1995, PASP, 107, 896 



-25- 



Shields, J. C, 1992, ApJ, 399, L27 

Shields, J. C, & Filippenko, A. V. 1990a, ApJ, 353, L7 

Shields, J. C, & Filippenko, A. V. 1990b, A J, 100, 1034 

Shopbell, P. L., & Bland-Hawthorn, J. 1998, ApJ, 493, 129 

ShuU, J. M., & McKee, C. 1979, ApJ, 227, 131 

Sparks, W. 1992, ApJ, 399, 66 

Storey, P. J., Mason, H. E., & Saraph, H. E., 1996, A&A, 309, 677 

Soker, N., & Sarazin, C. L. 1990, ApJ, 348, 73 

Taylor, G. B., Barton, E. J., & Ge, J. P. 1994, AJ, 107, 1942 

Unger, S. W., Taylor, K., Pedlar, A., Ghantaure, H. S., Penston, M. V., & Robinson, A. 
1990, MNRAS, 242, 33 

van Breugel, W., Heckman, T. M., & Miley, G. 1984, ApJ, 276, 79 

Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295 

Voit, G. M., & Donahue, M. 1997, ApJ, 486, 242 

Voit, G. M., & Donahue, M. 1990, ApJ, 360, L15 

Voit, G. M., & Donahue, M., & Slavin, J. D. 1994, ApJS, 95, 87 

Wirth, A., Kcnyon, S. J., & Hunter, D. A. 1983, ApJ, 269, 102 

Yan, L., & Cohen, J. G. 1995, ApJ, 454, 44 



This preprint was prepared with the AAS IMj7;X macros v4.0. 




Fig. 1. — Continuum-subtracted Hq;+[N II] image of NGC 1275 obtained on 1992 January 
2 U.T. by Shields and R. Pogge at the Perkins 1.8-m telescope, showing the position angles 
of the slit employed for the measurements presented here. 
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REST WAVELENGTH (A) 

Fig. 2. — Composite spectrum generated from averaging filament spectra obtained at PAs 
62°, 77°, and 283°. The emission features seen at 5057 A and 6629 A are [O III]A5007 and 
Ha, respectively, from the high-velocity system, redshifted by cz 3000 km s~^ relative to 
the low-velocity system. 
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Fig. 3. — Line-ratio diagrams displaying reddening-corrected values measured at different 
locations within the filaments. Approximate loci employed for classification of nebular types 
are indicated. Representative la error bars are shown. 
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Fig. 4. — Line-ratio diagram employing lines of oxygen only, with reddening-corrected 
filament measurements shown. The effect of reddening by Ay — 1 mag is given by the 



arrow. 
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Fig. 5. — Radial profiles of the [O IIl]/H/3 line ratio at the four PAs; error bars represent 
la uncertainties. 
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Fig. 6. — Expanded view of the composite spectrum from Figure 2, illustrating the mea- 
surement for NGC 1275 (solid histogram) and best-fitting continuum template (dashed line) 
in the vicinity of He 11 A4686. The expected location of this emission feature is shown. 
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Fig. 7. — Reddening-corrected [N II]A6583/Hq; line ratio as a function of Ha line width for 

the filament measurements. Also shown by the solid lines are predictions of shock models 
taken from Dopita & Sutherland (1995; DS = shock only, DSp = shock + precursor), Binette 
et al. (1985; BDT), and Shull & McKee (1979; SMK). In plotting the model predictions, we 
have equated the shock velocity to the FWHM of Ha. 
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Fig. 8. — Similar to Figure 6, illustrating the spectral region around [O III]A4363. 
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Fig. 9. — Model continua employed in photoionization calculations. The dotted line is the 
continuum incident on the screen, while the solid line is the transmitted continuum incident 
on the filament. Free-free, free-bound, and two-photon emission from the screen make the 
transmitted continuum higher than the incident one for energies below ~ 0.8 Ryd. 
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Table 1: Observed Emission-Line Fluxes 



H ,Tyi 1 GGi /"in 1 .in 






[U iiJAo/2/ 


z.5U±U.iU 






U. i4iztU.uo 


1 Q-UO 0/1 




n 90-1-0 0^ 


974-0 OA 


H7 


0.35±0.04 


0.42±0.05 


[0 III]A4363 


<0.035 


<0.04 


He II A4686 


<0.02 


<0.024 


H/3 


1.00±0.03 


1.00±0.03 


[0 III]A4959 


0.24±0.03 


0.23±0.03 


[0 III]A5007 


0.66±0.04 


0.62±0.04 


[N I]A5200 


0.33±0.03 


0.30±0.03 


He I A5876 


0.25±0.02 


0.19±0.02 


[0 I]A6300 


0.82±0.04 


0.59±0.03 


[0 I]A6364 


0.33±0.04 


0.24±0.03 


[N II]A6548 


1.22±0.06 


0.84±0.04 


Ha 


4.15±0.14 


2.86±0.10 


[N II]A6583 


3.54±0.12 


2.43±0.08 


He I A6678 


0.07±0.03 


0.05±0.02 


[S II]A6717 


1.38±0.06 


0.91±0.04 


[S II]A6731 


1.07±0.05 


0.70±0.03 



^ Normalized to H/3 = 1.93 x 10 erg s ^ cm ^ arcsec ^. 

^ Normalized to H/3 = 7.11 x 10~^^ erg cm~^ arcsec"^, after correction for Ay — 1.2 mag. 
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Table 2: Upper Limits on Coronal Lines and Theoretical Predictions 



Line/Ha^ 


Upper Limit^ 


Model^ 


Notes 


[Ni XII]A4232 


1.8x10-2 


1.2x10-^ 




[Ar XIV]A4414 


1.0x10-2 


1.0x10-3 


4 


[FeXIV]A5303 


9.4x10-3 


1.6x10-3 


5 


[Ca XV]A5445 


1.5x10-2 


0.6x10-3 




[Fe X]A6374 


1.4x10-2 


2.1x10-3 




[Ni XV]A6702 


7.9x10-3 


0.1x10-3 





^ Surface brightness of Ha is 1.9x10"^^ erg s"^ cm-2 arcsec-2^ after correction for Ay — 1.2 
mag. 

2 After correction for Ay = 1.2 mag. 

3 From Sarazin & Graney (1991), unless otherwise noted. 
^ From Voit et al. (1994). 

5 Voit et al. (1994) predict 15.6x10-3. 



Table 3: Line Ratios of the Photoionization+Screen Model 



Line Ratio 


Photoionization Model 


Observed Range 


He II A4686/H/3 


0.024 


<0.024 


He I A5876/H/3 


0.18 


0.16-0.20 


[0 l]A6300/H/3 


0.78 


0.33-0.78 


[N II]A6583/H/3 


2.59 


2.10-3.00 


[S II]A6717,6731/H/5 


3.00 


1.38-1.83 


[0 III]A5007/H/3 


0.43 


0.28-0.63 



